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Abstract. In this work we present a new Monte Carlo generator for Direct top and Single top 
production via flavour-changing neutral currents (FCNC). This new tool calculates the cross 
section and generates events with Next-to-Leading order precision for the Direct top process and 
Leading-Order precision for all other FCNC single top processes. A set of independent dimension 
six FCNC operators has been implemented - including four-fermion operators - where at least 
one top-quark is present in the interaction. 



1. Introduction 

With the new experimental data from the Large Hadron CoUider (LHC) one will be able to 
scrutinise the Standard Model (SM) boundaries with unprecedented precision. Flavour-changing 
neutral currents (FCNC) top decays such as t ^ V q, where V = Z,'y, g and q = u,c are highly 
suppressed in the SM due to the well known Glashow-Iliopoulos-Maiani (GIM) mechanism [1]. 
top-quark FCNC branching ratios are of the order of 10~^^ and 10^^^ for the electroweak and 
strong case respectively. However, some versions of two-Higgs doublet models (2HDMs) predict 
quite larger values. In fact, for some regions of the 2HDM parameter space the BR(t g q) 
could reach values of the order of 10~^ 0, which is eight orders of magnitude above the SM 
value. The high statistics expected in LHC data forces an increase in the precision of all relevant 
theoretical calculations. Hence, at least Next-to-Leading (NLO) calculations should be made 
available whenever possible. The experimental searches performed both at LHC and at the 
Tevatron on the direct top FCNC process have already used NLO cross sections from [3]. 
However, the NLO calculation was used only as a normalization factor and the events were 
generated using LO generators like TopRex ^ or Protos [5]. In this work we present a new 
Monte Carlo tool that generates FCNC direct top events with NLO precision. The FCNC 
interactions were implemented via an effective lagrangian [6] where a complete set of dimension 
six operators was used. The events are generated in the Les-Houches Event (LHE) format [7] 
and can therefore be easily interfaced with shower algorithms such as the one in Pythia [H]. 
Since the generator calculates the total inclusive cross section, the FCNC single top process was 
also included and can be generated independently. Not only strong FCNC interactions, but also 
electroweak and four-fermions dimensions six operators were included. Each operator can be 
turned on independently allowing for independent studies of the different Lorentz structures. In 



section [2] the effective formalism is presented while in section [3] we show the physical processes 
already included in MEtop. In section |4] we give a description of the algorithm used by the 
Monte Carlo generator in order to obtain a final NLO result. In section [5] we discuss a few 
selected results and in section [6] we sum up our conclusions. 

2. Effective operators 

In the effective lagrangian formalism we assume that the SM is no more than a low energy limit 
of some more general theory to be revealed at TeV scale. In this context, the top-quark FCNC 
interaction is seen at low energies as a contact interaction that can be parametrized with higher 
dimension effective operators. The complete lagrangian is written as 

A A2 

where we have truncated the expansion series at dimension d = 6. We follow the formalism in [6j, 
where the SM symmetries are imposed and all operators are built using the already known SM 
fields. This results in a reduced list of the possible dimension six operators. Additionally, 
baryon and lepton number conservation is imposed, which restricts the operators to dimension 
six only. Even imposing the SM symmetries, the number of all possible FCNC operators is 
quite vast. However, the list can be reduced by using appropriate equations of motion and 
Fierz transformations applied to the physical operators we are dealing with. The final minimal 
complete set of FCNC operators for top-quark physics was obtained in [HI [lOl EH 112] • With this 
procedure, a single operator survives in the strong sector [llj. Following the notation in [6j this 
operator can be written as 

0% = glA'^a^^n^^cAG"^^ (2) 

where G^^ is the gluonic field tensor, right-handed u or c-quark field and is the left-handed 
quark doublet. The same operator also appears in the literature as a dimension 5 operator 

iA^,|g-,AV^^(/. + /i.75)tG;i, , (3) 

where Kj is a real parameter, i = u,c, gs is the strong coupling constant and fi and hi are 
complex numbers which respect + = 1. This operator is responsible for the effective 
FCNC vertices tgu and tgc. The minimum set of operators in the electroweak sector is 

O'^B^ = Qi ^'^^ < 4> B,. , 0%,^ = qi n < 4> , (4) 



Ot = i{<l^'^D,<t>){uWu'n), (5) 
O^ll''^ = z {c^^D,^) (gl j^qi), Of/^ = i (0t rj D,^) {qi ^ tj gi) , (6) 

0%={<PU){qlu^^), (7) 

where B'^'^ and W^^^ are the U{1)y and SU{2)l field tensors, respectively. The field (/) is the 
usual Higgs boson doublet. In addition to this electroweak set of operators, we have included 
subsets of four-fermion operators which are described in [13] . 



3. Physical Processes 

The physical processes included in MEtop are FCNC direct top production and FCNC single 
top production. Direct top is available both at LO and at NLO while FCNC single top is availble 
only at LO. All processes can be generated for both pp and pp hadron colliders. Furthermore, the 
full top decay chain is included so that spin correlations are preserved. Also, in order to perform 
spin correlations studies, a on/off switch for the top decay chain was made available. In figures[T] 
and [2] we show the generic Feynman diagrams for the physical processes implemented in MEtop. 
The production process depicted in figure [T] has only strong FCNC operator contributions. Born 
level FCNC direct top is shown in figure [T] (left) together with the t+g process (right). In figure[2] 
we present the diagrams for the LO t + q production process from gluon fusion (right) and with 
quarks in the initial state (left) . The former is only affected by strong FCNC operators while 
the later has contributions from the strong, electroweak and four-fermion sectors. 
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Figure 1. FCNC leading order direct top production and top + gluon production at the parton 
level. Only FCNC strong operators contribute to the process. 
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Figure 2. FCNC leading order top + quark production at the parton level. FCNC strong and 
electroweak operators contribute to the process together with 4F operators. 



4. NLO approximation 

As mentioned in the previous section, MEtop is able to generate FCNC direct top events 
with NLO precision. The implementation was performed by adopting an NLO effective 
approximation [H]. A general NLO calculation of an n-particles final state process can be 
represented by 

da NLO = dasd^n + dayd^n + da^d^n+i, (8) 

where B stands for Born, V for Virtual and R for Real. It is well known that the virtual term 
has infrared divergences that will only be cured by the inclusion of the real radiation term aR. 
We show in figure [3] the class of diagrams contributing to NLO direct top production. The 
first process qg ^ t correspond to the Born term. The second, qg — >• t, corresponds to the 
Virtual term while the third, qg — >• gt, is the real radiation term. This last term gives rise 
to an infrared divergence which correspond to a low Pt and/or to a collinear external gluon. 
Kinematically, the divergence can be parametrized by the gluon transverse momentum, Pt, and 




Figure 3. Inclusive FCNC direct top production at NLO in QCD. 



it obviously occurs in the limit Pt — )• 0. In order to avoid the divergence one could just impose 
a Pt cut. The resulting region has to be excluded from the event space generation with some 
care - since the problem arises only in the collinear/soft limit, we can fill this phase space region 
with a branching mechanism. By doing so, we assume a collinear factorization, where the QCD 
radiation is emitted by one of the external legs of the Born configuration. This approximation 
can be translated by 

\M,+,\'d^2 ^ m'd^,^^P,,,,{z)dz^ (9) 

ZTT t ZTT 

where Mt+g and Mt are the real radiation for (t + g) and direct top amplitudes, respectively; d^i 
is the phase space for the ith body processes and Pq^qg{z) is the Altarelli-Parisi splitting function. 
The variable t is a resolution parameter which in our case is the top quark transverse momentum 
p!^^ . The hard Pt region will then be filled by using the appropriate transition amplitude and 
the low Pt region by the use of a shower mechanism. In this NLO effective approximation |14| . 
the Pt phase-space is then divided in two regions parametrized by which we call p^°'^'^^. 
This parameter plays the role of a matching parameter where the criterion for a good transition 
is the smoothness of the final Pt distribution. The final cross section free of divergences can be 
written as 

UNLO = KasiPP < Pr'"") + MPt > P^'^'^'^), (10) 

where Pj^^ is the transverse momentum of the parton shower emission. The contribution from 
the virtual term day in equation [8] is included via a K-factor applied to the born term and 
because born and virtual terms share the same 2 — )• 1 topology we assume they have similar 



kinematics. It is from equation 10 that events in MEtop are generated. Since MEtop is a parton 
level generator, the events must then be submitted to a parton shower simulator in order to 
include initial (ISR) and final (FSR) state radiation. Therefore, since these parton showers are 
based on a branching mechanism, we leave the collinear and/or soft Pt phase space region to 
be filled by the corresponding parton shower. In order to avoid double counting we just have to 
ensure that the parton shower will not populate the Pt > Pi^"'^'^^ region in the first emission. 
To prevent it, the shower mechanism must be performed with the PT-ordered scheme with the 
first emission starting at Pt = For the NLO direct top process, MEtop assumes that 

the Pythia shower will be used with its PT-ordered scheme [15] flag on. Finally, we note that 
because this is a 2 — )• 1 process, the shower mechanism implemented in Pythia will not add FSR 
but only ISR. This does not pose any problems because FSR is negligible when compared to 
ISR due the large top-quark mass. 



5. The generation process and some results 

In figure 4 we show the diagram flow of the MEtop generation process. It was drawn specificahy 
for the NLO direct top case. The remaining processes follow the same flow except for the 
K-factor calculation step. The amplitudes used to evaluate the cross sections were generated 
with CalcHEP |16j . whit Feynman rules derived with LanHEP [17J and all integrations were 
performed using the Cuba Library [18]. The K-factors are calculated "on the fly" and therefore 
each sub process has its own normalization factor. 
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Figure 4. MEtop flow cart for NLO Direct top generation. 



In figure [5] we present the top-quark Pt distribution for direct top after the first emission 
and with the starting scale shower at mt (black solid line). The blue dashed line represents 
the top Pt distribution of the real radiation process. As explained before, both regions overlap 
and therefore a matching must be introduced in order to avoid double counting. The matching 
process must be chosen so that the transition between the rescaled direct top distribution and the 
real radiation process is smooth. In figure|6]we show the top Pt distribution for p^"-^'^^ = IQ, 15 
Gev where direct top is the gray dashed line and t + is the blue dashed line. The first one 
corresponds to the top distribution after the first ISR emission, where the PT-shower starting 
scale was set to 10 and 15 Gev, respectively; the second one is the top distribution from the 't-|-g' 
process where a P^"* of 10 and 15 Gev was imposed. We have studied several p™"*^'* scenarios 
in the 5 Gev to 20 Gev range finding no significant differences. We therefore recommend that 
this matching value should be included in any analysis as a systematic error. 

In figure [T] we show the final NLO direct top quark Pt and rj distribution (solid line), as well 
as the LO result (dashed line) after full ISR, FSR and Multiple Interactions (MP). It is clear 
that in both distributions the final NLO result does not amount to simply multiplying the LO 
result by a K-factor. In fact, the distributions show that the NLO contribution concentrates 
more events in the low Pt region and produces top quarks at higher angles. In figure [8] we 
show the equivalent distribution for the lepton coming from the decaying W. The 6-quark and 
neutrino distributions show similar differences between the NLO and LO results. For the top 
decay products, the NLO contribution does not result in a major change in distributions making 
a K-factor adjustment suitable. However, one must pay special attention to analysis where the 
top quark momentum is reconstructed. For example, analysis which include the study of the 
top spin polarisation usually requires full reconstruction of its momentum. In this case, figure [7] 
show that NLO events should be considered. Finally, when studying the inclusive NLO direct 
top production one must add the t + jet (single top FCNC) process as well. The distribution of 
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Figure 5. Pt distribution of the top quark for ^/s = 7 TeV. The black soUd hne is for direct 
top production after the first branching in ISR, with starting scale of mj. The blue dashed line 
is for the hard process top+gluon production. 
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Figure 6. Pt distribution of top quark after the first ISR branching with a p^°-^'^^ of 10 GeV 
(left) and 15 GeV (right). 
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Figure 7. Comparison of the LO and NLO Pt (left) and r] (right) distributions of the top 
quark at the partonic level after the full shower (ISR+FSR) and Multiple Interaction. 



the sum of the two processes is shown in figure [9] for the top Pt (right) and the top (left). The 
FCNC Single top process was generated with a Pj!^^ = 10 Gev and only subprocesses where a 
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Figure 8. Comparison of the LO and NLO Pt (left) and i] (right) distributions of the lepton 
from t — )• bW — )• hlv at the partonic level after the full shower (ISR+FSR) and Multiple 
Interaction. 



FCNC interaction takes place were included. Processes like ud — )■ th which are pure SM single 
top production were discarded. In tables [T] and [2] we show the total contributions from each 
LO and NLO processes. They were calculated for ku^c = 0.01 and for LHC @ 7 TeV. As we 
can see from figure [9j the FCNC single top contribution does not produce considerable shape 
modifications, although its contribution to the total cross section is approximately 13%. 
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Figure 9. Pt (left) and r] (right) distributions of the top quark for NLO direct top (solid line) 
and NLO direct top plus pp tq with ^ iq Q^y ^^^^ > iq Q^y 



Table 1. FCNC Direct top (anti-top) LO and NLO cross sections for ku,c = 0.01 and LHC @ 
7 Tev. 



Subprocess 


LO 


NLO 


ug ^ t 


6.12 


8.74 


eg t 


0.91 


1.67 



6. Conclusions 

We have presented a new Monte Carlo generator dedicated to top FCNC physics. We have 
included the main FCNC top production channels: direct top and single top. The direct top 



Table 2. FCNC t + q with ku,c = 0.01 and LHC @ 7 Tev and PTcut = 10 Gev. 



Subprocess ku = 0.01 kc = 0.01 
p,p^t,q 1.12 0.40 



channel is available at NLO and the single top at LO. A complete set of dimension six operators 
are already included. We have shown that the NLO top distributions obtained for the direct 
top cannot be built via a K-factor from the LO ones. Therefore NLO events should be used 
especially in the cases where top-quark momentum reconstruction is required. 
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